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The three-dimensional momenta of quarks inside a hadron are encoded in transverse momentum- 
dependent parton distribution functions (TMDs). This work presents an exploratory lattice QCD 
study of a TMD observable in the pion describing the Boer-Mulders effect, which is related to po¬ 
larized quark transverse momentum in an unpolarized hadron. Particular emphasis is placed on the 
behavior as a function of a Collins-Soper evolution parameter quantifying the relative rapidity of 
the struck quark and the initial hadron, e.g., in a semi-inclusive deep inelastic scattering (SIDIS) 
process. The lattice calculation, performed at the pion mass m,, = 518 MeV, utilizes a definition of 
TMDs via hadronic matrix elements of a quark bilocal operator with a staple-shaped gauge connec¬ 
tion; in this context, the evolution parameter is related to the staple direction. By parametrizing 
the aforementioned matrix elements in terms of invariant amplitudes, the problem can be cast in a 
Lorentz frame suited for the lattice calculation. In contrast to an earlier nucleon study, due to the 
lower mass of the pion, the calculated data enable quantitative statements about the physically in¬ 
teresting limit of large relative rapidity. In passing, the similarity between the Boer-Mulders effects 
extracted in the pion and the nucleon is noted. 


I. INTRODUCTION 

Transverse momentum-dependent parton distributions (TMDs) [1] constitute one of the pillars on which the three- 
dimensional tomography of hadrons rests. Together with the three-dimensional spatial information derived from 
generalized parton distributions (GPDs), they permit a comprehensive reconstruction of hadron substructure and 
thus have a bearing on seminal topics in hadron physics, such as orbital angular momentum contributions to nucleon 
spin, or spin-orbit correlations in hadrons. Through the selection of particular parton spin and transverse momentum 
components, a variety of effects can be probed, including naively time-reversal odd (T-odd) quantities such as the 
Sivers and Boer-Mulders functions; these only exist by virtue of initial or final state interactions in corresponding 
physical processes, introducing a preferred chronology in the description of the process. For example, in semi-inclusive 
deep inelastic scattering (SIDIS), the operative element is final-state interactions between the struck quark and the 
hadron remnant; on the other hand, in the Drell-Yan (DY) process, initial state interactions before the lepton pair 
production enable T-odd effects. TMDs thus in general have to be considered in the context of specific physical 
processes, within a factorization framework appropriate for the process in question, separating the hard reaction from 
the TMD and other elements such as fragmentation functions. In the case of T-odd effects, the process-dependence 
manifests itself in the prediction of a sign change of the Sivers and Boer-Mulders functions between the SIDIS and 
DY processes [2]. 

In view of the fundamental importance of TMDs and the rich spectrum of effects that can be probed, TMDs 
have been, and continue to be the target of a variety of experimental efforts. Deep-inelastic scattering experiments 
performed by COMPASS [3], HERMES [1] and Jefferson Lab [5] have yielded TMD data including evidence for the 
T-odd Sivers effect. Complementary Drell-Yan experiments at COMPASS [5] and Fermilab [7] are envisaged, which 
could, in particular, test the aforementioned sign change between the SIDIS and DY processes. Related transverse 
single-spin asymmetries have been measured at RHIC in polarized proton-proton collisions [S] . Further experimental 
efforts at RHIC are projected to provide insight into strong QCD evolution effects expected for the Sivers TMD [5]. 
TMDs furthermore constitute a central focus of the proposed Electron-Ion Collider facility m- 

To complement these efforts, providing nonperturbative QCD input from first principles to the analysis of TMD 
effects, a project to calculate TMD observables within lattice QCD was initiated and developed in [TTHT3] . The 
present work constitutes a continuation of this project. As described in detail below, the formal definition of TMDs 
is based on nonlocal operators, specifically quark bilocal operators with a gauge connection that takes the shape 
of a staple. The path followed by the gauge connection is in principle infinite in length, and thus it cannot be 
straightforwardly treated in terms of an operator product expansion, as is commonly done, e.g., for ordinary parton 
distributions or generalized parton distributions. In view of this genuinely nonlocal character of the operators, lattice 
QCD explorations of corresponding hadronic matrix elements directly at the nonlocal operator level were undertaken 
in [mill], concentrating initially on the simpler case of straight gauge links connecting the quark operators. The 
nonlocal nature of the operators in particular raises novel questions regarding regularization and renormalization, 
which were addressed in considerable detail in [12j . Whereas these questions deserve further study, the aforementioned 
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explorations suggest that it is a viable working assumption to treat nonlocal lattice operators in analogy to the fashion 
in which they are treated in continuum QCD m, namely, by absorbing divergences into multiplicative soft factors. 
These soft factors can then be canceled in appropriate ratios; this scheme was used to construct TMD observables in 
the subsequent investigation m, and will be used in the present work. Formally related studies of nonlocal lattice 
operators, in which a gauge link in the (Euclidean) time direction originates from the propagation of a heavy auxiliary 
quark, have been carried out in [15]; also, a direct approach to light-cone distribution amplitudes based on nonlocal 
lattice operators was laid out in m- Moreover, the comprehensive framework for investigating parton physics within 
lattice QCD put forward in m and developed and explored in [TBH^ relies on a direct treatment of such nonlocal 
lattice operators. 

The present work focuses on a TMD observable related to the Boer-Mulders effect in a pion. Lattice QCD studies 
of pion structure, predominantly focusing on form factors, have been previously reported in Choosing the 

pion as the hadron state is motivated by the principal goal of the investigation presented here, namely, understanding 
the behavior of TMDs as a function of an evolution parameter quantifying the rapidity difference between the hadron 
momentum and a vector describing the trajectory of the struck quark. Details are furnished further below. In the 
previous nucleon study |13] . no definite conclusions regarding the limit of large rapidity difference proved possible. By 
virtue of its lower mass, the pion provides a larger rapidity difference at given momentum, and this choice of hadron 
state thus aids in approaching the limit of physical interest. In addition, the spinless nature of the pion permits 
additional spatial averaging to suppress statistical uncertainties. Indeed, the chief advance of the present work lies 
in providing quantitative insight into the limit of large evolution parameter. Preliminary accounts of this work were 
given in [5UH5^ . 


II. DEFINITION OF TMD OBSERVABLES 
A. Correlation functions 

Quark transverse momentum-dependent parton distributions (TMDs) can be defined in terms of the fundamental 
correlator 


^Ssubtr.(^-P>---) = \{P\ 9(0) TU[0,r]v,T]v + b,b] q{b) \P) 


( 1 ) 


where P denotes the momentum of the hadron state; the present work focuses on pions, and thus no spin is attached to 
the state. F represents an arbitrary Dirac 7 matrix structure. The quark operators at positions 0 and b are connected 
by the gauge link U[0,r]v,r]v + 6 , 6 ], which connects the points listed in its argument by straight-line segments; thus, 
the gauge link has the shape of a staple, cf. Fig. with the unit vector v specifying the staple direction and t] its 
length. One is ultimately interested in the limit 77 —>■ 00 , which in a concrete lattice calculation is of course reached 
by extrapolation. This gauge link form incorporates final state interactions between the struck quark and the hadron 
remnant in semi-inclusive deep inelastic scattering (SIDIS) [33] . and analogously initial state interactions in the Drell- 
Yan process (DY). The ellipsis in the argument of ‘hJujgujjtr indicates that the correlator will depend on various further 
parameters, related, e.g., to regularization, specified below as needed. 

Fourier transformation of Q, 


^^^\x,kT,P, ... 
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FIG. 1: Gauge link structure 6/[0, rjv, rjv + 6 , 6 ] in the correlator 0 - In a concrete lattice calculation, the limit 7 —>■ 00 has to 
be taken numerically. 
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leads to the momentum space correlator which ultimately will be parametrized in terms of TMDs, cf. below. 
The position space correlator Q, as written, requires regularization not only of the quark operator self-energies, but 
also of the self-energy of the Wilson line (this is indicated by the subscript “unsubtr.”). The regularization of the 
latter is effected by dividing by the soft factor S. The detailed structure of the soft factor depends on the concrete 
factorization approach employed. For example, in the scheme developed in HUSH, which, for reasons discussed 
further below, provides the phenomenological framework for the present study, the soft factor takes the form 






S{0)ibT,+oo,-oo)S(^0){bT,ys,-oo) 

S(o){bT,+oo,ys) 


(3) 


where 5 '(o)(&T) ■ • •) is a vacuum expectation value of Wilson line structures extending, initially, into space-like direc¬ 
tions; some of them remain at finite rapidity j/g, whereas others are taken to the light-cone limit, i.e., infinite rapidity. 
This particular form of the soft factor, containing more than one rapidity, cannot be cast in a Lorentz frame in 
which it exists at a single time, and for this reason, it is not suited for evaluation within lattice QCD. However, the 
observables that will be defined further below are ratios in which the soft factors cancel [13]. There is, therefore, no 
obstacle to evaluating those observables within lattice QCD, and the detailed form of the soft factor is immaterial. 
An alternative construction of TMD soft factors which, in principle, is amenable to lattice QCD evaluation has been 
put forward in j 191 within the framework laid out in m- 

As written in the transverse components br of the quark separation b are Fourier conjugate to the quark 
transverse momentum kx, while the longitudinal component b - P is Fourier conjugate to the longitudinal momentum 
fraction x. The present work will be confined to the case b ■ P = 0, corresponding to evaluating the integral with 
respect to x over the correlator and TMDs derived from it. It should be stressed, however, that there is no 
obstacle to extending calculations of the type presented here to a scan of the b ■ P-dependence^, yielding upon Fourier 
transformation the cc-dependence of $1^1 and the TMDs under consideration. Studies of the b ■ P-dependence in the 
straight gauge link {ijv = 0) case have already been carried out in |12j , and further investigations in this direction are 
planned for future work. A related proposal to obtain the cc-dependence of parton distributions has been put forward 
and explored in mHio]. 

Finally, ([^ is evaluated at b'^ = 0, in accordance with the standard phenomenological framework, which employs 
a Lorentz frame in which the hadron of mass nih propagates with a large momentum in 3-direction, P'^ = (P*^ -|- 
P^)j\^ ^ ruh] then, the quark momentum components scale such that the correlator ([^ and TMDs derived from 
it are principally functions of the quark longitudinal momentum fraction x = k^/P'^ and the quark transverse 
momentum vector kr, whereas the dependence on the component k~ = {k^ — k'^)/\/2 <C becomes irrelevant in 
this limit. Correspondingly, ([^ is regarded as having been integrated over implying that the conjugate variable 
6+ is to be set to zero, as written. 

Before continuing with the parametrization of in terms of TMDs, it is important to note a further specification 
with respect to the staple direction v. The Wilson lines along the legs of the staple represent an effective, resummed 
description of gluon exchanges between, in the case of SIDIS, the struck quark and the hadron remnant [33]. Accord¬ 
ingly, their direction v should be taken to follow the path of the ejected quark, close to the light cone from the point of 
view of the hadron. Whereas, at tree level, there is no obstacle to the most straightforward choice, namely, a light-like 
V, beyond tree level, this choice is associated with rapidity divergences [3S]. Various schemes have been advanced to 
treat these divergences [IllllllEelET], and the equivalence between some of them discussed in [38] . In particular, the 
scheme advanced in [HIM] effects regularization by tilting the staple direction v slightly off the light cone into the 
spacelike region. This feature is crucial for a concrete implementation of a lattice QCD evaluation of matrix elements 
of the type Q, as will be detailed further below. Thus, the phenomenological framework providing the backdrop 
for the present treatment is, specifically, the one advanced in [HIM]- Choosing v to be spacelike for the purposes 
of regularization implies a dependence of the calculation on an additional parameter, characterizing proximity of the 
staple to the light cone. Here, this parameter will be chosen as 


c 


V ■ P 

RM ’ 


(4) 


where the absolute length of a four-vector w is denoted by \w\ = \/|wR. In terms of C, the light-cone limit corresponds 


^ In a practical calculation, the range of accessible b- P is limited by the available b and P, \b‘ P\ < \P\V —6^ (where P denotes the spatial 
momentum), leading to an increasing systematic uncertainty at small x. 
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to C —>■ oo. The generated lattice data thus have to be extrapolated to the double limit ry —>■ oo, C —)■ oo, i.e., the limit 
of infinite staple length, and staple direction converging toward the light cone. 

Alternative to this purely kinematic characterization, also the parameter ^ = 2mhC is frequently employed [39], 
and viewed as a dynamical scale, to be compared to Aqcd] perturbative evolution equations in ( can then be derived 
at sufficiently large cf., e.g., [ID]. The question whether one has reached the asymptotic regime appropriate for 
the definition of TMDs presumably has both kinematic and dynamical aspects. The applicability of perturbation 
theory, e.g., for determining evolution equations is a dynamical issue most adequately characterized by considering 
the dimensionful parameter in relation to the characteristic QCD scale, i.e., requiring ^ ^ ^qcd- On the other 
hand, this by itself does not necessarily guarantee kinematics close to the light cone, C ^ 1- Presumably, both 
conditions need to be taken into account in general. In the present work, no quantitative connection to perturbative 
evolution is attempted; the numerical data do not reach values of ^ which lie clearly within the perturbative regime. 
The maximum spatial momentum employed is 1.17 GeV, and the maximum value of C is C = 2.1 GeV, corresponding to 
C = 2.03. Instead, the dependence of the lattice data on the kinematical variable ^ will be studied empirically, including 
exploration of ad hoc ansatze for the large-C behavior, allowing for corresponding extrapolations. Gharacterizing the 
large-C limit is in fact the primary goal of the present investigation. This limit was seen to present a considerable 
challenge in the previous study m, which considered nucleon TMDs; no definite statements concerning the large-^ 


behavior proved possible. The present work, focusing on pions, permits accessing higher C both by virtue of the 
lighter hadron mass (note that the hadron mass enters the denominator of and by employing additional spatial 
averaging facilitated by the spinless nature of the pion, enhancing statistics. As will be seen further below, the data 
extracted in the present investigation are of sufficient quality to yield a signal for the C —>■ oo limit of the generalized 
Boer-Mulders shift defined in eq. (23) below. This constitutes the main advance of this work. 


B. Parametrizations 


Returning to the momentum space correlator 
yields, at leading twist. 


its parametrization in terms of the relevant Lorentz structures 




(5) 

( 6 ) 


For spinless particles such as the pion, there are only two leading twist TMDs, in contrast to the eight which arise 
for spin-1 particles |lT1143j . The TMD /i is simply the unpolarized quark distribution, whereas the Boer-Mulders 
function [33] encodes the distribution of transversely polarized quarks in the pion. The Boer-Mulders function is 
odd under time reversal (T-odd). Physically, it only exists by virtue of the final and initial state interactions, in the 
SIDIS and DY processes, respectively, which break the symmetry of the processes under time reversal. Formally, it 
is the introduction of the additional vector v describing the staple direction in the staple-shaped gauge link which 
breaks the symmetry; the Boer-Mulders function vanishes for a straight gauge link, rjv = 0. Correspondingly, the 
generalized Boer-Mulders shift defined in eq. (23) below will be an odd function of v. 

On the other hand, one can also decompose the position space correlator into invariant amplitudes [13] . 

The general decomposition is (the combinations corresponding specifically to the leading twist TMDs ([^ and ([^ will 
be considered further below): 


2 unsubtr. 

II 

(7) 

2 unsubtr. 

2 

'ryi ^ ^ 

= P^A2 imllPA3+ "v^Bi 

V ■ F 

(8) 

2 ^unsubtr. 

2 

7 777 

= P^bpV^Bi 

V ‘ F 

(9) 

2 ^unsubtr. 

Q 

^ TYl —’ ITTl'^ 

V • F V ‘ F 

(10) 


The present treatment focuses on the special case b ■ P = 0, which in the context of TMDs, defined in a frame in 
which 6+ = 0 and vt = Pt = 0, also implies b ■ v = 0 m- Under these constraints, the above relations are readily 
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inverted. The amplitudes needed below are, explicitly, 


^2 

A 4 
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( 11 ) 

( 12 ) 

(13) 

(14) 


Note that Bi and B 3 are regular for v ■ P —>■ 0 owing to the prefactor. Of particular interest are the leading twist 
objects 


1 

1 

2 ]P^ 


$ 


[7+] 

unsubtr. 


$ 


unsubtr. 


= A 2 B 
= —irriT^e 


(15) 

(16) 


(where 7,j denote transverse spatial indices), given in terms of the amplitude combinations 


A 2 B — A 2 + ^1 — 1 + C Bi 

Aab = ^4 — ^1 — 1 + 


(17) 

(18) 


Also these combinations are regular for C —)■ 0 by virtue of the prefactors in (12),(14). Note that, in the case 
of vanishing spatial momentum, one cannot identify “forward” and “backward” directions for u; there is then only 
a single br anc h in \rjv\, the sign of which is a matter of definition. Although only the combinations A 2 B and Aab 
appear in (15),(16), for the numerical analysis to follow, it will be valuable to be able to consider A 2 , Bi, Aa and B 3 
individually, not just those combinations. 

Given that and are related via a Fourier transformation, the quantities arising in the respective 

decompositions ([^-([^ and (15)-(16) must be similarly related, i.e., the amplitudes AiB must be related to Fourier- 
transformed TMDs. Indeed, denoting x-moments of generic Fourier-transformed TMDs by 


/HG)(6^, ...)=„! dxx^-^ J d^kr /(x, ...) (19) 

= (1^) M\bT\\kT\)f{x,k^,...) (20) 

where J„ denotes the Bessel functions of the first kind, one finds [13] 

fP°\b^T,t...,Vv-P) = 2A2B{-b^T.b-P = 0,b-v = 0,C,7jvP)/S{b^T,...) (21) 

/i^W«(6^,C,...,r?u-P) = 2AAB{-b^T,b-P = 0,b-v = 0,C,7jvP)/S{b^T,...) (22) 


Note the appearance of the soft factors on the right hand sides. 


C. Boer-Mulders shift 


As already indicated further above, one obtains an observable in which the soft factors cancel by forming a suitable 
ratio, namely, the “generalized Boer-Mulders shift” 


{ky)uT{brp, ...) = 


__ hp^^\b'^,...) _ A4b(-5|,0,0,C,7?u ■ P) 

A2b(- 6|,0, 0,C,r7v • P) 


/P^(fe^,...) 


= m. 


(23) 
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Note that ratios of this type also cancel F-independent multiplicative field renormalization constants attached to the 
quark operators in Q at finite physical separation b. It should be emphasized that the construction thus far is a 
continuum perturbative QCD construction. That this construction carries across into lattice QCD, i.e., that the 
lattice operators are similarly regularized and renormalized by multiplicative soft factors which cancel in ratios, is a 
working assumption which was already explored in considerable detail in [12) . and which will be investigated further 
in future work. Physically, this assumption appears plausible at least at separations substantially larger than the 
lattice spacing, where the lattice operators are expected to approximate the corresponding continuum operators. 

To interpret the generalized Boer-Mulders shift, note that the br —>■ 0 limit of the quantities defined in (20) formally 
corresponds to fc^-moments of TMDs, 


/M(")(o,...) = I'^dxx^-^ J d^kr 


Thus, in the formal 

{ky)uT{Q, ■ ■ ■) = 


0 limit, the generalized Boer-Mulders shift reduces to the “Boer-Mulders shift” 
^^[^'^^^(0,...) _ f dx f d^kr kr, P,...) 


7™^°^ (0, . . .) Jdxf d?kT + (a;, kT,P,...) 


(24) 


(25) 


ST = (1,0) 


which, in view of the structure of the right-hand side, formally takes the form of the average transverse momentum 
in y-direction of quarks polarized in the transverse (“T”) cc-direction, in an unpolarized (“17”) pion, normalized to 
the corresponding number of valence quarks. This is in accord with the (s^ x ba’) • P structure of the correlator 
(16), where st represents the quark spin. The numerator in (25) sums over the contributions from quarks and 


antiquarks, whereas the denominator contains the difference between quark and antiquark contribu tion s, thus giving 
the number of valence quarks [I^ I45j . It should be noted, however, that the fc7"™oments of TMDs (24) appearing in 
(25) are in general divergent |46j at large kx and thus not well-defined absent further regularization. The generalized 


quantity (23) is a natural regularization, with finite bx effectively acting as a regulator through the associated Bessel 


weighting, cf. (20). This Bessel weighting also is advantageous in the analysis of experimental asymmetries 


In the present work, lattice QCD data for the generalized Boer-Mulders shift (23) will be obtained and presented at 


finite br- The path by which these data can be obtained proceeds via lattice QCD evaluation of the fundamental 
correlator ([^ for a range of Dirac and staple link structures, extraction of the relevant invariant amplitudes ©-(HI, 
and construction of the ratio ( [2^ . As already mentioned further above, for this lattice QCD calculational scheme to 
be viable, it is necessary to employ a phenomenological framework such as the one advanced in [HIM], in which all 
separations in the correlator ([^ are spacelike, including the staple direction v. Such a scheme implies dependence 
on the Collins-Soper-type evolution parameter Q cf. (j^, quantifying proximity of the staple to the light cone. The 
principal focus of the present investigation is, indeed, the dependence of the generalized Boer-Mulders shift on 
including its asymptotic ^—>"00 behavior. 


III. LATTICE QCD CALCULATIONS 

Lattice QCD employs a Euclidean time coordinate, serving to project out hadronic ground states via the associated 
exponentially decaying time evolution. As a consequence, when evaluating matrix elements of operators in hadronic 
states, no Minkowski time separations in those operators can be accomodated; one is restricted to operators which 
are defined at one single time. This is the reason why it is imperative to employ a framework in which all separations 
in the fundamental correlator ([^ are spacelike. Only in this case is there no obstacle to boosting the problem to 
a Lorentz frame in which the operator in 0 exists at a single time, and performing the lattice calculation in that 
particular frame. 

The decomposition of the resulting correlators into invariant amplitudes, cf. 0 -(jldj), is a further crucial element of 
the present treatment. Expressed in this fashion, the results of the lattice calculation are immediately applicable also 
in the original Lorentz frame in which ([^ was initially defined. Einally, as already discussed above, the construction 
of ratios of amplitudes in which soft factors cancel serves to connect the results to phenomenological observables such 
as the generalized Boer-Mulders shift ( |^ . 

The lattice QCD data for the present exploration were obtained within a mixed action scheme employing domain 
wall valence quarks on an Nf = 2 -f 1 dynamical asqtad quark gauge ensemble provided by the MILC collaboration 
[49| . Since the principal focus lies on understanding the systematics of the large limit, which proved inaccessible 
in previous investigations, a fairly high pion mass, rriT^ = 518 MeV, was chosen for this study to alleviate statistical 
fluctuations. Further details of the ensemble are given in Table |l| This mixed action scheme, including the specific 
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xT 

a(fm) 

amu,d 

arris 

(MeV) 

(GeV) 

#conf. 

#meas. 

20® X 64 

0.11849(14)(99) 

0.02 

0.05 

0 

06 

1.348(09)(13) 

486 

3888 


TABLE I: Parameters of the lattice ensemble. Note that the lattice spacing a was determined in a different scheme than 
in [SSI m], and consequently it, as well as the pion and nucleon masses quoted, differ slightly from the ones given in the 
aforementioned references m- The first error quoted for the pion and nucleon masses is statistical, the second stems from the 
conversion to physical units using a. The bare asqtad quark masses are denoted mu,d,s- Eight measurements were made on 
each gauge configuration. 


ensemble employed here, has been used extensively by the LHP Collaboration for studies of hadron structure, cf., 
e.g., [5U1IM] . It also provided the basis for the previous nucleon TMD investigation reported in m- 

To extract the correlator Q, one evaluates both three-point functions and two-point functions C-^pt with pion 
sources and sinks of definite spatial momentum^ P, 

X,- ,Xf 

X,- .X/ 

where U, t and tf are source time, operator insertion time and sink time, respectively, and (j) denotes an interpolating 
field with the quantum numbers of the pion. Wuppertal-smeared quark fields were employed to construct these pion 
sources and sinks, which were separated hy tf — = 9a. Only connected contractions contributing to C^pt were 

evaluated; disconnected contributions, which are expected to be small, are omitted in all results presented below. In 
the case of the tt'*' meson, w-quark and d-quark distributions coincide; contrary to the nucleon case, there is therefore 
no nontrivial u — d quark combination in which disconnected contributions exactly cancel. 

The fundamental correlator Q is then obtained from plateaus in t for t tf in the three-point to two-point 
function ratio m, 


(26) 

(27) 


*[r] _ p(p\ ^3pt [O] {ti,t,tf,P) 

^ ^ C2p,{U,tf,P) 


(28) 


where E{P) is the energy of the pion state and O is taken to be the operator in ([^. It should be noted that, at the 
employed source-sink separation oi tf — ti = 9a = 1.07 fm, signihcant excited state contaminations in the plateaus 
extracted from ( |28[ ) cannot be excluded a priori. This issue was not investigated in the present exploratory study 
at the fairly high pion mass = 518 MeV. However, in future work at lower pion masses, where excited state 
contaminations are exacerbated, it will present an additional challenge. 

The set of combinations of pion momenta and staple-shaped gauge link paths used is listed in Table [H] It should 
be noted that, in the case of either b or v extending into a direction which does not coincide with a lattice axis, 
there is more than one optimal approximation of the corresponding continuum path by a lattice link path; e.g., if 
b = 2 (ei -I- 62), where e.^ denotes the lattice link vector in ^direction, both the sequence of links (01,02,61,62) and 

the sequence of links ( 02 , 61 , 62 , 61 ) equally well approximate the continuum path. In such a case, $(i„subtr always 
averaged over all equivalent lattice link paths, for both b and v vectors. This symmetry improvement of the lattice 
operators is important to preserve the manifest time-reversal transformation properties present for the continuum 
staple-shaped gauge link path operators. 

Note also that, in the mixed action scheme used for these calculations, before evaluating domain wall propagators 
for valence quarks, the asqtad gauge configurations are HYP-smeared to reduce dislocations (or rough fields) that 
would otherwise allow right-handed states on one domain wall to mix with left-handed states on the other domain 
wall. The lattice gauge link paths in Q were constructed using those same HYP-smeared gauge configurations. This 
has the advantageous consequence that renormalization constants are closer to their tree-level values, while it would 
have no effect in the continuum limit. On the other hand, as shown in ref. |20j . there are significant differences 
between 0 and 2 steps of HYP-smearing in the direct calculation of parton distributions using straight gauge link 
paths instead of staples (see Figs. 3-5 therein), so the optimal use of HYP smearing requires further study. 


^ In practice, momentum conservation eliminates the need for projection at the source, provided one projects onto zero momentum transfer 
at the operator insertion in Cspt instead. 













h/a 

rj\! a 

P . aLj{2'K) 

n ■ (0,0,1), n = —9,..., 9 

in' • (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0), 

(0,-2,0), (-1,-1,0), (-2,-1,0) 

in' ■ (0,1,0) 

(0,0,0), (-1,0,0), (-2,0,0), 

(0,-2,0), (-1,-1,0), (-2,-1,0) 

in'.(1,1,0) 

(0,0,0), (-1,0,0), (-2,0,0), 

(0,-2,0), (-1,-1,0), (-2,-1,0) 

in'.(1,-1,0) 

(0,0,0), (-1,0,0), (-2,0,0), (0,-2,0), (-2,-l,0) 

in'.(2,1,0) 

(-2,-l,0) 

in'.(2,-1,0) 

(-2,-l,0) 

n ■ (0,1, 0), n = —9,..., 9 

in' ■ (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0), (0,0,-2) 

in'. (0,0,1) 

(0,0,0), (-1,0,0), (-2,0,0), (0,0,-2) 

in'. (1,0,1) 

(0,0,0), (-1,0,0), (-2,0,0), (0,0,-2) 

in'. (1,0,-1) 

(0,0,0), (-1,0,0), (-2,0,0), (0,0,-2) 

n ■ (1,0, 0), n = —9,..., 9 

in'. (0,1,0) 

(0,-2,0), (0,0,-2) 

in'. (0,0,1) 

(0,-2,0), (0,0,-2) 

in'. (0,1,1) 

(0,-2,0), (0,0,-2) 

in'. (0,1,-1) 

(0,-2,0), (0,0,-2) 

n- (0,1, l),n = -4,... ,4 

in' . (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0) 

n- (0, l,-l),n = -4,... ,4 

in' . (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0), (-1,-1,-!) 

in'. (0,1,1) 

(-1,-1,-!) 

n- (1,0, l),n = -4,... ,4 

in' . (0,1,0) 

(0,-2,0) 

n- (1,0,-l),n = -4,... ,4 

in'. (0,1,0) 

(0,-2,0), (-1,-1,-!) 

in'. (1,0,1) 

(-1,-1,-!) 

n- (l,l,0),n = -4,... ,4 

in'. (0,0,1) 

(0,0,-2) 

n- (l,-l,0),n = -4,... ,4 

in'. (0,0,1) 

(0,0,-2), (-1,-1,0), (-1,-1,-!) 

in'. (1,1,0) 

(-1,-1,0), (-1,-1,-!) 

n ■ (0, 2,1), n = —3,..., 3 

in' . (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0) 

n • (0, 2, — 1), n = —3,..., 3 

in' . (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0) 

n ■ (0,1, 2), n = —3,..., 3 

in' . (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0) 

n ■ (0,1, —2), n = —3,..., 3 

in' . (1,0,0) 

(0,0,0), (-1,0,0), (-2,0,0) 

n ■ (2,0,1), n = —3,..., 3 

in'. (0,1,0) 

(0,-2,0) 

n • (2,0, — 1), n = —3,..., 3 

in' . (0,1,0) 

(0,-2,0) 

n ■ (1,0, 2), n = —3,..., 3 

in'. (0,1,0) 

(0,-2,0) 

n ■ (1,0, —2), n = —3,..., 3 

in' . (0, 1,0) 

(0,-2,0) 

n ■ (2,1, 0), n = —3,..., 3 

in'. (0,0,1) 

(0,0,-2) 

n • (2, — 1,0), n = —3,..., 3 

in'. (0,0,1) 

(0,0,-2) 

n ■ (1, 2, 0), n = —3,..., 3 

in'. (0,0,1) 

(0,0,-2) 

n • (1, —2,0), n = —3,..., 3 

in'. (0,0,1) 

(0,0,-2), (-2,-1,0) 

in'. (2,1,0) 

(-2,-l,0) 


TABLE II: Sets of staple-shaped gauge link paths and spatial pion momenta P used on the lattice. Gauge link paths are 
characterized by the quark separation vector b and the staple vector jyv, cf. Fig. The surveyed range of r], parameterized 
in the table by the integer n', was always chosen to extend from zero to well beyond the point where a numerical signal 
ceases to be discernible. The maximal magnitude of the Collins-Soper parameter ( attained in these sets is |(^| = 2.03, for 
P • aL/ {2n) = (—2,0, 0) paired with rivja = in' • (1, 0,0) and rotations thereof. 
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FIG. 2: Generalized Boer-Mulders shift for u-quarks as a function of staple extent ri\v\, for fixed ^ = 1.01; the panels illustrate 
data obtained at a succession of quark separations |feT|. Plateau fits and extraction of asymptotic values (open symbols) are 
described in the main text. 


IV. NUMERICAL RESULTS 
A. SIDIS and DY limits 

The first step in the analysis of the obtained data concerns the behavior as a function of staple length r]. For 
ease of notation, both positive and negative rj are considered for a fixed u • P > 0 to distinguish staples oriented 
in the forward and backward directions with respect to the pion momentum. Of particular physical interest is the 
asymptotic behavior for t] ±oo, corresponding to the SIDIS and DY limits. Fig. [^displays results for the u-quark 
generalized Boer-Mulders shift as a function of 77|z;| at a fixed ^ = 1.01, with each of the four panels corresponding to 
a successively larger transverse quark separation {BtI- The T-odd behavior of the observable is evident. As the SIDIS 
and DY limits are approached, a clear plateau behavior in r;|n| is observed up to moderate values of \bT\; as |6 t| rises, 
statistical uncertainties increase (as indicated by the jackknife error estimates in the plots), and the identification of 
the plateaus becomes more tenuous, cf., e.g., the lower right panel, for |6t| = 0.48 fm. Plateau values are extracted 
by averaging over the regions 7a < r]\v\ < 9a and —7a > r]\v\ > —9a, respectively, as indicated by the fit lines in 
the plots; finally, the SIDIS and DY limits are obtained imposing T-oddness, i.e., the two plateau values in each plot 
are averaged with a relative minus sign to yield the asymptotic SIDIS and DY estimates also displayed in the panels 
(open symbols). The asymptotic values slightly decrease in magnitude as |6 t| rises. 

Fig. summarizes the results obtained in the SIDIS limit as a function of the quark separation Ibx], for three 
different values of the Collins-Soper parameter (. Note that the data at small |6 t|, up to |6 t| « 2a = 0.24 fm, may be 
affected by discretization artefacts, but at larger |6j’|, the data are expected to well approximate the continuum limit. 
For larger cf. also further examples below, the statistical fluctuations rapidly increase, and no useful signal was 
obtained beyond |6 t| = 0.5 fm in the case ( = 2.03. The data appear to approach well-defined limits as either \bT\ or ( 
becomes large. The behavior as a function of C will be discussed in greater detail below; the behavior as |5 t| becomes 
large seems plausible: Physically, once |6 t| exceeds the size of the pion, the correlator Q cannot anymore probe 


















10 


> 

lU 

o 


H 



0.0 

-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1— 

--'-1-^ 


Boer-Mulders Shift (SIDIS), u-quarks 

-O.I 

T 


- 0.2 

: ijii It 





-0.3 

• 



^ ■ 

^=2.03 : 

-0.4 

• ♦ 

^= 1.01 : 


, = 518 Me V • 

^=0 

-0.5 

' . ■ 1 ■ ■ . 1 ■ ■ . r , , , 1 , 

^_._I_^ 


5 ^ ' 

0.0 0.2 0.4 0.6 0.8 1.0 


Ibrl (fm) 


FIG. 3: Generalized Boer-Mulders shift for u-quarks in the SIDIS limit as a function of the quark separation {brl, for three 
different values of the Gollins-Soper parameter The shown uncertainties are statistical jackknife errors. 


correlations inside the pion; it only contains vacuum-vacuum and vacuum-pion correlations. The 1 67 ’!-dependence 
of these correlations is then expected to be dominated by the typical exponential fall-off with |6t| observed in the 
vacuum^; neglecting all other dependences in comparison, and canceling the dominant behavior in the ratio (23) leads 
to the expectation of a constant asymptotic behavior in I^tI- 

As C increases, one furthermore would expect the aforementioned constant to converge to an asymptotic limit. 
Remarkably, however, even the C = 0 data appear to already approach the same large-| 6 t| constant as the data 
at higher Possibly, this may be understood as a consequence of the Lorentz invariance of the vacuum; only 
vacuum-vacuum and vacuum-pion correlations are probed at large l^-rl, and it seems plausible that these would be 
independent of the pion momentum. It would be desirable to gain a more definite understanding of this property. 
The most remarkable feature of the data, however, is the apparent tendency of the generalized Boer-Mulders shift to 
become constant in |6t| as C is increased, not only for asymptotic values of |6t|, but for all |6t|- No obvious reason for 
this behavior at low to intermediate {bxl is apparent, and it would be very interesting to develop an understanding of 
it. The constant behavior implies that, in the relevant transverse momentum |A:t| range corresponding to the probed 
range of |&t|, the transverse momentum spectrum of polarized quarks is the same as for unpolarized ones. 


B. Evolution in ^ 

A special focus of the present investigation is the behavior as a function of C and the large C limit, i.e., studying 
in detail the behavior of the sequence of data seen in Fig. [^at a fixed value of | 6 t|- Fig. shows two further panels 
analogous to the ones in Fig. but with ( varying between the panels and {BtI = 0.36 fm fixed instead. Thus, in 
terms of a ^ sequence, the lower left panel of Fig. 0 lies in between the two panels displayed in Fig. |4| As already 
mentioned further above, for pion spatial momentum P = 0, corre^onding to the left panel in Fig. I^there is only 
one branch as a function of ri\v\, as shown. The right panel in Fig. W corresponding to C = 2.03, illustrates the rapid 
deterioration of signal as the pion momentum is increased. Nevertheless, at the moderate | 6 t| used here, a plateau 
can still be extracted. 

Fig. El summarizes all results obtained in the SIDIS limit at |6t| = 0.36 fm as a function of the Collins-Soper 
parameter ( (left panel), and analogous data for the nearby value |6t| = 0.34fm (right panel). From the figure, the 


® Heuristically, the expectation value of the gauge link staple, which, after integrating out the quark fields, may be thought of as being 
completed into a closed loop by the (fluctuating) world lines of dynamically propagating quarks, is expected to be determined by the 
chromodynamic flux piercing the loop. This is, e.g., the origin of the Wilson loop area law demonstrating conflnement in Yang-Mills 
theory m- Once exceeds the size of the pion, at most one of the legs of the staple can traverse the pion; the other runs entirely 
within the vacuum. Consider now varying by shifting the latter leg; the area being added or removed from the loop lies purely 
within the vacuum. Only vacuum chromodynamic flux is being added or subtracted, while the chromodynamic flux influenced by the 
pion remains flxed. Thus, the variation of the expectation value with |6x| is determined purely by vacuum properties (which is not to say 
that the expectation value becomes entirely independent of the properties of the pion; only its |6(r|-dependence does). This argument 
is unchanged if one subsequently averages over different positions of the leg of the staple traversing the pion. 
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FIG. 5: Generalized Boer-Mulders shift for u-quarks in the SIDIS limit as a function of the Gollins-Soper parameter for 
fixed quark separation |6t| = 0.36fm (left) and |6t| = 0.34fm (right). Filled symbols correspond to the full generalized 
Boer-Mulders shift, whereas open symbols correspond to the partial contribution obtained by replacing A 4 B in ( |23| by A 4 , for 
further discussion, cf. main text. Spatial pion momenta P along various directions, as labeled, were employed; labels in general 
subsume more than one magnitude along the specified axis, cf. Table [11] for the complete set of momenta for which data were 
generated. Data labeled P ~ (1,0,0) include averaging over 90 degree rotations thereof, i.e., momenta along any of the three 
spatial lattice axes are included. A given value of ^ can be accessed using different pairs of P and staple direction v; in such 
cases, data corresponding to different P directions are slightly offset from one another horizontally for better visibility. The 
shown uncertainties are statistical jackknife errors. 


good rotational properties of the calculation are evident; a given value of C can be accessed using different directions 
of the pion spatial momentum P and the staple direction v, including both on- and off-axis directions, as shown. The 
corresponding results coincide, indicating that potential lattice artefacts are under control, and thus buttressing the 
physical significance of the data obtained. To assess the asymptotic behavior at large (, it is advantageous to consider 
not only the full generalized Boer-Mulders shift, but also the partial contribution obtained W replacing A 4 B in (23) 
by A4, omitting the contribution from i? 3 , cf. eq. (18). Both quantities are displayed Fig. as is evident from the 
figure, the partial A 4 contribution vanishes at (^ = 0, but monotonically increases in magnitude as ( is increased. By 
contrast, the remaining contribution from B 3 to the full generalized Boer-Mulders shift, dominant at (^ = 0, decreases 
in magnitude as ( rises. This matches the behavior expected from eq. (18), according to which the contribution from 
B 3 becomes insignificant for (( —>■ 00 under the assumption that the amplitude B 3 stays hnite. Altogether, the full 
generalized Boer-Mulders shift also decreases in magnitude as ( increases. Thus, the behavior of the data suggests 
that, by considering both quantities, one has access to both lower and upper bounds for the generalized Boer-Mulders 
shift, considerably increasing the confidence in the extrapolations to large ( discussed below. The comparison between 
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FIG. 6: Least-squares fits to the generalized Boer-Mulders shift data shown in Fig. (left). Separate fits to filled and open 
data points are shown. Left panel displays fits obtained with ansatz c -|- d/(; right panel displays fits obtained with ansatz 
c-I- d/G- Error bands show statistical jackknife uncertainties of the fits. Data points that lie outside the horizontal range of 
the displayed fit curves were not included in fits; in the case of the open symbols, one less data point was taken into account 
because it already lies in the region of an apparent inflection point, at which the simple asymptotic fit ansatz certainly is not 
applicable anymore. The asymptotic values displayed include the separate fits to the filled and open data points (filled and 
open squares), as well as the results of combined fits (filled circles); for details, cf. main text. The shown uncertainties are 
statistical jackknife errors. 
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FIG. 7: Least-squares fits to the generalized Boer-Mulders shift data shown in Fig. (right). Separate fits to filled and open 
data points are shown. Left panel displays fits obtained with ansatz c -I- d/(; right panel displays fits obtained with ansatz 
c-I- d/G- Error bands show statistical jackknife uncertainties of the fits. Data points that lie outside the horizontal range of 
the displayed fit curves were not included in fits. The asymptotic values displayed include the separate fits to the filled and 
open data points (filled and open squares), as well as the results of combined fits (filled circles); for details, cf. main text. The 
shown uncertainties are statistical jackknife errors. 


the partial and full quantities also permits an assessment of the extent to which evolution in ( has progressed towards 
the asymptotic limit. Evidently, according to Fig. already about half of the magnitude of the full generalized 
Boer-Mulders shift is subsumed in the partial A 4 contribution at C « 2. A significant part of the ( evolution has thus 
already been achieved at that value of <^. In this respect, the present study yields a much clearer picture than was 
obtained in the previous nucleon investigation |13] . 

To obtain quantitative statements about the large ( limit, least-squares fits^ to the ^-dependences of the data were 


4 


Whereas the data at higher generally display the larger statistical uncertainties, the data at lower ^ are likely to deviate systematically 
to a larger degree from any putative simple asymptotic behavior. It therefore does not seem appropriate to bias the fits towards the 
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FIG. 8: Least-squares fits to generalized Boer-Mulders shift data at |6t| = 0.27fm, analogous to Figs.|^and[^ The data were 
all obtained using spatial pion momenta P along the lattice axes (P ~ (1, 0, 0) in the nomenclature of Fig.j^ including multiples 
and 90 degree rotations thereof, cf. Table [11] for the complete set of momenta for which data were generated). Separate fits to 
filled and open data points are shown. Left panel displays fits obtained with ansatz c + d/(; right panel displays fits obtained 
with ansatz c + d/(^. Error bands show statistical jackknife uncertainties of the fits. Data points that lie outside the horizontal 
range of the displayed fit curves were not included in fits. The asymptotic values displayed include the separate fits to the filled 
and open data points (filled and open squares), as well as the results of combined fits (filled circles); for details, cf. main text. 
The shown uncertainties are statistical jackknife errors. 



Fit function 

Full B.-M. shift 
(GeV) 

Gontribution A4 
only (GeV) 

Combined fit 
(GeV) 

RMS deviation of 
combined fit (GeV) 

|6t| = 0.36 fm 

c -1- d/( 

-0.146(26) 

-0.141(36) 

-0.145(25) 

0.00755 

|6t| = 0.36 fm 

c + d/C^ 

-0.166(16) 

-0.110(22) 

-0.148(15) 

0.01695 

|&t| = 0.34 fm 

c + d/C 

-0.145(33) 

-0.112(33) 

-0.128(29) 

0.01466 

|6t| = 0.34 fm 

c + d/C^ 

-0.157(19) 

-0.084(19) 

-0.121(16) 

0.02315 

|6t| = 0.27 fm 

c + d/C 

-0.153(37) 

-0.122(29) 

-0.138(28) 

0.00975 

|6t| = 0.27 fm 

c + d/C^ 

-0.170(25) 

-0.095(20) 

-0.133(19) 

0.03855 


TABLE III: Asymptotic values of least-squares fits to the ^-dependences of the full generalized Boer-Mulders shift and the 
partial A 4 contribution separately, as well as a combined fit. These are the values displayed at 1 ^ = 00 in Eigs.jbl^ Uncertainties 
quoted are statistical jackknife errors. The right-hand column gives the root mean square deviation, per degree of freedom, of 
the combined fit from the data means. 


performed, using power-correct ion ansatze of the form c -I- d/( and c -I- d/C^. It should be noted that the values of ( 
for which data were obtained in the present study do not reach clearly into the perturbative regime within which a 
reliable connection to perturbative evolution can be guaranteed. The aforementioned fit functions should therefore be 
regarded as ad hoc ansatze. On the one hand, separate fits to the full generalized Boer-Mulders shift data and to the 
partial A 4 contribution were performed; on the other hand, a combined fit to both quantities with a common constant 
c (but, of course, separate coefficients d) was carried out. The results are shown in Figs. 6][8 and the asymptotic 
values summarized in Table [Tn| Fi g. J6| displays the fits to the data of Fig. (left), at |6 t| = 0.36 fm, whereas Fig. [7] 
displays the fits to the data of Fig. [5](right), at |6t| = 0.34fm. Fig. [^additionally displays fits to data obtained at a 


lower |6t|, namely, |6t| = 0.27fm. 

The fits using the form c + d/C are superior to those using the form c ■ 


d/C^. In the former case, the asymptotic 


values obtained by analyzing the full generalized Boer-Mulders shift data and the partial A 4 contribution separately 
agree within uncertainties, and also with the result of the combined fit. In the latter case, the asymptotic values 
obtained with the separate fits differ significantly. Not surprisingly, the results of the combined fits using the two 
fit ansatze are quite similar; by construction, they settle roughly in the middle between the two separate quantities. 


lower C data by performing a fit taking into account the statistical uncertainties. 
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Nevertheless, also in the combined fits, the c + d/C ansatz is more favorable, as evidenced by the last column in 
Table III which gives the root mean square deviation, per degree of freedom, of the combined fit from the data 
means. The deviation is considerably larger for the c + d/C^ fits; it is of the order of, or larger than, the jackknife 


statistical uncertainty of the asymptotic value. By contrast, for the c + d/C ansatz, the deviation is considerably 
smaller than the jackknife statistical uncertainty of the asymptotic value. These observations favor the c + d/C fit 
ansatz. 

Note that the asymptotic values for the generalized Boer-Mulders shift obtained using the c + d/C ansatz both at 
{brl = 0.36fm and at {brl = 0.34fm, as well as at the lower separation \bT\ = 0.27fm, all coincide within uncertainties. 
These extrapolations thus do not modify the observation made further above in connection with Fig. namely, that 
the generalized Boer-Mulders shift appears to become constant in I&tI as C becomes large. 

Altogether, it is encouraging to observe that the data generated in the present study are of sufficient quality to 
ensure that the signal survives extrapolation to the C oo limit, with an uncertainty as low as, roughly, 20% for 
the selected parameter values analyzed above. It is possible to tentatively discriminate between different asymptotic 
models. Further improvements of the analysis are to be expected by confronting lattice data with perturbative 
evolution equations for the observables considered here. 


C. Pion-proton comparison 

In [^, the question was posed how similar Boer-Mulders functions are in different hadrons, and also for different 
flavors. The present study permits addressing this question in passing, since Boer-Mulders data are now available 
on the same lattice ensemble both for nucleons [13] and for pions, obtained here, with corresponding parameters, 
facilitating a comparison. Fig. provides a sample juxtaposition of the tt“'"- meson and proton u-quark generalized 
Boer-Mulders shifts generated on the same ensemble, at identical transverse quark separation [bxl = 0.36fm, and 
identical spatial hadron momentum P and staple direction v, as a function of staple extent. This corresponds to 
identical C = but of course C then differs because of the appearance of the hadron mass in the denominator 

of Q. In this particular juxtaposition, the behavior of the two shifts is seen to be very similar. In the proton, the 
approach to the large \r]\ plateaus appears to be slightly faster; the plateaus lie within 10% of one another. One could 
argue that the more appropriate comparison is one at identical and not at identical Cj as shown here. However, 
using identical makes the correspondence hardly less favorable. If one interpolates the SIDIS limit pion data at 
Ibxl = 0.36fm shown in Fig.(left) from ( = 1.01 to ^ = 0.39, the ( value corresponding to the proton data in Fig. 
(right), then the pion data are enhanced in magnitude by around 20%, rendering the pion data about 10% higher in 
magnitude than the proton data, as opposed to 10% lower in the comparison at equal 

Of course, it should be noted that a close quantitative correspondence to the degree observed here for the u-quark 
generalized Boer-Mulders shifts is not generic. In the tt “'“-meson, the d-quark Boer-Mulders shift is identical to the u- 
quark one, whereas in the proton, the d-quark Boer-Mulders shift is appreciably higher in magnitude compared to the 
M-quark one m- The corresponding comparison would therefore not reveal a similarly close quantitative agreement. 
Nevertheless, qualitative features such as the signs of the different flavor contributions remain in correspondence. 


V. SUMMARY AND OUTLOOK 


Building on the previous study [13] of transverse momentum-dependent parton distributions in the nucleon, the 
present work focused in particular on the evolution of TMD ratios such as the generalized Boer-Mulders shift (23) 
as a function of the Collins-Soper parameter cf. (|^. To this end, the behavior of this TMD observable in a pion 
was explored; the pion, by virtue of its lower mass, allows one to access higher and the spinless nature of the pion 
facilitates spatial averaging to increase statistics. Indeed, in contrast to the nucleon study [13], which did not yield 
definite conclusions concerning the large-/' limit, the present investigation provides data of sufficient quality to perform 
tentative extrapolations to / = oo. While the generated data do not reach into the / region in which a clear connection 
to perturbative evolution equations is given, the ad hoc extrapolations explored here lead to estimates of the ( = oo 
limit of the generalized Boer-Mulders shift ( |2^ to within an uncertainty of as low as roughly 20%, for favorable values 
of the transverse quark separation br, cf. Table m The confidence in these extrapolations is buttressed in particular 
by the ability to partition the generalized Boer-Mulders shift into two separate contributions, as discussed in detail 
in section IV B[ such that one is provided with both an upper and a lower bound for the asymptotic behavior. 

Having obtained data for the generalized Boer-Mulders shift in the pion, the present study, together with [T3|, in 
passing also permits a juxtaposition of pion and nucleon TMD observables. A close correspondence between the u- 
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FIG. 9: Generalized Boer-Mulders shift for u-quarks in a tt"'"- meson (left) and a proton |13| (right), as a function of staple 
extent ri\v\, for quark separation |6 t| = 0.36fm. The data in the two panels were obtained using the same spatial hadron 
momentum P and staple direction v. This corresponds to identical ( — 2mhC, but, consequently, differing (, as labeled in the 
panels, cf. main text for further discussion. Plateau fits and extraction of asymptotic values are analogous to Fig. In the 
proton case, taken from the analysis described in m, a wider fit range was used. 


> 0-4 

o 

D 

0.2 

5 0.0 
- 0.2 
S -0.4 


<S> 


Boer-Mulders Shift, u-quarks 


(= 1.01 
Ibj-I = 0.36 fm 
ffi;, = 518MeV 

E— DY 




SIDIS 


- 

Boer-Mulders Shift 

- 

^ 1 




♦ 

up - quarks, * -r 


- 

^ = 0.39, ’ • • 

z ’ 


|b| = 0.36fm 



DY SIDIS —» 


— 00 

-10 -50 5 10 

oo 


Staple extent n|v| (lattice units) 


quark generalized Boer-Mulders shifts in a Tr'^-meson and a proton is observed, corroborating the behavior conjectured 

in [Ml[53]. 

Going forward, the success in determining the characteristics of the ^ evolution of the generalized Boer-Mulders 
shift in a pion suggests that analogous results will also be accessible for nucleons in future lattice studies with higher 
statistics, focused on obtaining signals at higher nucleon momenta. To this end, it will be useful to explore schemes 
of generating nucleon interpolating fields with a favorable overlap with higher momentum states. Further efforts 
currently in progress are concerned with the universality of TMD ratios of the type considered here under variation of 
the lattice discretization scheme; this is expected to provide empirical support for the working assumption employed 
here, that the lattice TMD operators are regularized and renormalized in analogy to the corresponding continuum 
operators, by multiplicative soft factors which cancel in TMD ratios. Moreover, lattice TMD calculations are in 
progress at lighter pion masses, with a view to transferring the exploratory results obtained here to the physical light 
quark mass regime. 
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